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•  ttaanlccad  Th«  factor  DT  «tap«ada  oa  nha  Incidant  wBYofroot, 

eha  curvacurao  of  tha  latarfacot  and  cha  mlacla%.  Indleaa  of  eha  t«o  nadla. 
Wo  axplielt  natrla  fonulaa  for  caleulaclns  llluatraca  lea  phyalcal 
aifniflcanca  via  ananplaa,  and  point  out  an  arronaoua  aolutlon  In  tha 
meant  litaratura. 

Mast,  wa  atu^tha  traaaniaaion  of  a  apharlcal  alactronafnatlc  wava 
threttfh  a  dialactric  ohall.  Tha  cuo  aurfaeaa  of  tha  ahall  ara  apharlcal 
(althar  coneava  or  convas),  and  chair  eantara  ara  arbitrarily  locatad  In 
mlatlon  to  Cha  aourea  point.  Tha  flald  aolutlon  datamlnad  by  cha 
gaonatrical  optica  thaory  la  glvan  In  a  alaipla  eloaad  fora.  Spaclal 
attantlon  la  glvan  to  cha  lana  efface  of  tha  dialactric  ohall  which 
convatta  tha  inconlng  apharlcal  pancll  Into  a  focualng  pancll. 

Finally,  wa  praaant  axtanalva  nuaarlcal  raaulta  on  arraya  o.*  point 
aourcaa  covarad  by  practical  radoaaa. 


□□ 


/J8TRACT 


In  •  prnvioua  report  (S],  dataila  of  Lha  analyala  of  a  curved  dlalactrlc 
radom  ualng  gaowatrlcal  optica  ware  conslaarad  for  a  point  aourca.  In  thia 
rapoit,  wa  firat  conaldar  the  tranaalaalon  of  a  apharical  or  plana  wave 
through  a  curved  dlalactrlc  Interface.  The  tranaalttad  field  la  proportional 
to  i.ha  product  of  the  conventional  Fraanal'a  transnlaalon  coefficient  and  a 
dl\«rgan(<*  factor  DF,  which  daacrlbaa  the  croaa-aactlonal  variation 
(convarganra  or  divergence)  of  a  ray  pencil  aa  the  latter  propagatea  In 
the  trananltted  region.  The  factor  DF  dependa  on  the  incident  wavefront, 
the  curvaturea  cf  the  Interface,  and  the  relative  Indicaa  of  the  two  aedla. 

We  give  explicit  aatrlx  fomitlaa  for  calculating  DF,  illustrate  Its  physical 
significance  via  exanples,  and  point  out  an  erroneous  nolutlrn  la  the 
recent  literature. 

Next,  we  study  Che  tranniiaslon  of  a  spherical  elactronagnetlc  wave 
through  a  dielectric  shell.  The  two  surfaces  of  the  shell  are  spherical 
(either  concave  or  convex) ,  and  their  centers  are  arbitrarily  located  In 
relation  to  the  source  point.  Tl'.e  field  solution  detemlned  by  the 
geoBctrlcal  optics  theory  Is  given  In  a  simple  closed  form.  Special 
attention  is  given  to  the  lens  effect  of  the  dielectric  shell  which 
converts  the  inconlng  spherical  pencil  into  a  focusing  pencil. 

FinalJv,  we  present  exterslve  nuserlcal  results  on  arrays  of  point 
sources  covered  by  practical  radones. 
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I .  IMTHOD'JCTION 


Many  practical  antannaa  ara  covarad  by  raUoaas,  whoaa  af facta  on 
tha  antanna  radiation  ara  of  conaidarabla  iaportanca,  aapacially  in 
today'*  hi^li-parforaanca  radar /coaminicat ion  ayataaa.  In  tha  paat 
quartar  of  a  cantury,  aavaral  atandard  analyaaa  hava  baan  daviaed  for 
analysing  radoca  af facta.  Nona  of  thaa  ia  axact,  and  laprovaaanta  ara 
alwaya  naadad.  Tha  praaant  raport  daacribad  an  affort  in  thia  diraction. 

A  typical  radoaa  problaa  aay  ba  atatad  aa  followa.  Lat  an  apartura 
antanna  A,  for  inatanca  a  horn,  a  alot,  or  a  conformal  array,  radiata  a 
known  fiald  {^(r)  in  fraa  apaca  (aaa  Flgura  ]).  A  protactiva  ahlald  or 
radoma  I  ia  placad  around  antanna  A.  Tha  proolam  ia  to  datarmina  tha 
radiation  fiald  f  for  tha  coapoaita  atructura,  l.a.,  tha  antanna  A 
radiating  in  tha  praaanca  of  tha  radoma.  Thia  problam  haa  racaivad  a  graat 
daal  of  attantion  from  many  raaaarchara  during  tha  laat  tvo  dacadaa,  and  a 
ao-callad  "bast  avallabla"  mathod  for  attacking  this  problam  appaara  to 
hava  amargad.  A  briaf  daacrlptlon  of  thia  mathod  ia  givan  b^low 

(a)  In  tha  vicinity  of  £,  tha  Incidant  fiald  ia  not  a  ray  fiald 
(locally  plana  wava) .  To  clrc*i*vant  thia  difficulty,  lat  bi  raaolv^id 
into  a  apactrum  of  plana  waves,  namely, 


i^r) 


.r 


Ik.  •  r 

dk_  1  dk^  W(k) 


a.i) 


'S(k)  ■  dx  j  dy  E^(r) 

Here,  Ic  «  (k  ,k  ,k  )  is  the  direction  of  propagation  of  the  piano-wave 
X  y  a 

spectral  component.  The  spectral  wave  number  in  the  z-direction. 


(1.2) 
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(a)  free  apace 


(b)  radome  environment 


Figure  1.  Antenna  A  and  radome  £. 


Figure  2.  Two  choices  of  Incident  directions:  A1  and  PI 


k  -  -  (k^  +  k^) 

*  *  y 


(1.3) 


may  be  real  (honogenous  plane  wave)  or  Imaginary  (inhomogenoua  plane 
wave) .  ‘the  weighting  factor  i((k)  is  the  amplitude  of  the  plane<>wave 
spectral  component  propagating  In  the  direction 

(b)  For  each  plane-wave  component,  a  tranomlsslon  coefficient  matrix 
TQ(k)  for  a  flat  dielectric  slab  can  be  obtained  from  any  staudard  text  on 
EM  theory.  The  subscript  zero  of  Indicates  that  it  Is  derived  from  the 
assumption  of  a  plane-wave  Incident  field.  The  transmitted  field  {^(?)  at 
point  2  on  the  outer  surface  of  Z  Is  calculated  from  the  formula 


1*^(2)  -  I  dk^  I*  dky  TqCU)  5(vt) 


(1.4) 


(c>  Once  £^(2)  la  known  for  all  points  on  the  outer  surface  of 
equivalent  surfoce  current  sources  (1(2),  ^(2))  can  be  determined.  The 
convolution  of  the  source  with  the  Green's  function  gives  the  desired 
radiation  field  which  is  expressible  as 


f 

t(r)  -  i  (G,  1  +  G,  it)  da 

}  I  1  ^ 


(1.5) 


outer  £ 


The  approach  described  above  Is  of  course  theoretically  sound.  However, 
its  faithful  execution  la  Impractical  becavise  of  the  extremely  labox'ious 
numerical  integrations  in  (1.4)  and  (1.5).  In  the  well-quoted  analyses 
by  Faria  [1]  and  Wu  and  Rudduck  [2],  the  numerical  Integration  In  (1.4) 
Is  avoided  hy  approximating  the  transmitted  field  at  point  2  by 


Ik  •  r 

S'=(2)  :  fQ({Q)  S(ltp)  e  ° 


(1.6a) 


where  the  Incident  dirnctlon  Is  determined  by 
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■  actual  ray  dlractlon  Al,  or  the  dlractlon  PI  of  tha 

Poynclng  vector  of  (Fifura  2).  (1.6b) 

Note  that  the  approximation  In  (1.6)  la  to  deacrlbe  by  a  plane  wave. 

Since  the  radome  la  In  the  near  xono  of  the  antenna,  thla  plane  wave 
approximation  for  deacrlbed  In  [1],  [2]  doea  not  seem  to  be  a  good  one. 

In  the  present  report,  we  approach  the  radome  problem  from  a  different 
viewpoint.  Instead  of  decompoalug  Che  incident  field  into  a  plane  wave 
spectrum,  we  approximate  the  flnlce-slzed  antenna  A  In  Figure  3a  by  an 
array  B  In  Figure  3b.  Each  element  In  array  B  radiates  a  spherical  wave. 
Those  spherical  wave  constituents,  transmitting  through  the  radome  £,  are 
superimposed  to  give  rlae  to  the  desired  radlaclou  field  ^  In  the  far  sone, 
Thus,  Che  key  step  in  Che  present  approach  Is  to  determine  the  transmission 
of  a  spherical  wave  through  a  curved  dielectric  shell. 

We  shall  apply  geometrical  optics  to  solve  the  transmission  through 
the  curved  radome.  Specifically,  we  first  study  in  detail  the  different 
aspects  of  refraction  at  a.  single  curved  dielectric  Interface;  second,  we 
analyze  the  behavior  of  a  spherical  dielectric  shell;  and  third,  we  apply 
our  theory  tc  various  practical  radomes. 
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Figure  3.  An  aperture  antenna  A  Inalde  a  radoM  la  approximated 
by  an  array  B.  Each  point  source  In  array  B  radiates 
a  spherical  wave. 


Figure  4.  Transmission  through  a  dielectric  shell  due  to  Incidence 
from  a  point  source  at  P^. 


II.  DBSCRIPTKMe  OF  PIIOBLBM 


7h*  MOMtry  of  th«  rodoMt  problMi  undar  conMldaratlon  la  akatchad 
In  Figure  A  point  aourea  at  producaa  a  apharlnal  vava  vhich  goaa 
through  a  curved  dielectric  ahall  with  nonunlfora  thlcknaaa.  Ray 
tachnlquaa  are  uaad  to  datarmlna  the  field  at  point  on  a  given  aurfaca 
outside  the  ahall.  First,  let  us  daacrlta  the  various  alaaents  Involved 
In  the  problem. 

Coordinate  Svatama  and  Tima  Convention,  the  auiln  coordinate  syatan 
la  the  rectangular  aystan  (x.y.s),  whose  origin  lo  chotan  at  the  source 
point  and  the  z>coordlnate  is  In  .a  direction  of  the  bean  maximum 
of  the  antenna.  Other  coordinate  aystama  at  points  F^,  P2  and  P^  along 
the  ray  are  defined  later.  The  field  la  tl»e-hanM>nle  with  the  time 
factor  expC-fJut)  which  la  suppreseed  throughout. 

Source .  Ve  essuma  Chet  the  source  has  a  well-defined  "phase  center" 

AAA 

at  point  P^,  the  origin  of  the  coordinate  system  (x.y.z),  and  radiates 
a  spherical  wave  denoted  by  (f^,  If  the  antenna  la  an  array  of  point 

sources,  it  is  necessary  to  consider  each  element  In  the  array  separately 
and  superlmpoae  their  final  fields  at  the  observation  points. 

Dlelnctric  radosw.  The  radome  la  a  dielectric  shell  with  nonuniform 
thickness  of  relative  dielectric  constant  •  c/c^  or  refraction  Index 
n  ■  1'^,  and  Is  bounded  by  the  inner  and  outer  surfaces  and 
respectively.  The  inner  surface  (near  the  source)  Is  described  by 
the  equation: 

z  -  fj^(x,y)  ,  for  Sj^  <  X  <  bj^  and  Cj^  <  y  <  dj^  .  (2.1) 

The  outer  surface  ^2  Riven  by  the  equation: 
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s  -  fjU.y) 


for  *2  '  ^2  ^2  ^  ^  ^  ^2 


(2.2) 


» 


It  is  not  nocossar/  to  Itnow  tho  wsnlytlcnl  fox«  of  tho  functions  fj^(x,y) 

and  f2(xty).  In  eonputstlon,  only  s  sot  of  discrsts  dsts  points 

[x  ,y  tf  }  vlch  n  «  1,2, ....M  is  nsodod  for  tho  doocrlptlon  of  f(f.  or  f,). 
n  n  n  x  « 

Thoso  points  sro  fitted  by  o  cub^ic  spllns  which  glvos  sutonstleolly  first 

snd  second  psrclsl  d^rlvstlvos  of  f,  l.o.,  9f/)s,  2f/2y,  9^f/9s^,  2^f/2x2y, 
2  2 

snd  3  f/3y  .  Thors  sro  two  ro^ulrononts  for  tho  cublc»splino  fltt 
(1)  tho  dots  points  esn  ho  dlstrlbutod  owor  s  rondos  t^ld,  but  thoy  oust 
bo  donso  onough  to  doscribo  thii  flno  dotslls  of  or  £2)* 
donsin  of  tho  dots  points  (s  <  x  <  b  snd  c  <  y  <  d)  wist  bo  sonawhst 
groator  than  tho  srss  of  £  In  which  tho  Incldoet  ray  Is  sxpoctsd  to 
Intorsoct  tho  rsdono. 

Obsorvstlon  points.  Obsorvotlon  point  lo  locstod  on  a  prospoclflod 
surfaco  £j,  which  can  bo  olthor  oao  of  tho  following  two  typos: 

(1)  Sphorlcal  £^  trltb  cantor  at  and  an  Inflnltoly  largo  radius. 

In  this  caso,  is  In  tho  far  flold,  and  tho  flold  at  P^ 
calculatsd  by  tho  ray  technlquo  Is  tho  final  rosulr . 

(11)  Planar  which  Is  just  outsido  tho  radono  and  nomol  to  tho 

c-axls.  In  this  caso,  wo  have  to  Intograta  tho  flold  on  to 
obtain  tho  fat  flold. 

In  lator  calculations,  wo  uso  itostly  tho  sphorlcal  £^  In  (1). 
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Ill.  FIELD  S^OLOTlOll  BT  CIONITEICAL  OPTICS 


For  •  giv«n  iocldooc  f'JLtld  t*^**'*^^  by  the  •ourct  at  point 

Pq  (Fig.  4).  Chn  •■ynptotlc  nolutlon  of  th«  finld  at  point  la  datar- 
mitiad  ualng  gaoMtrlral  optica  [3]t  (4).  The  Mthnd  of  solution  la 
daacrlbad  balow. 

A.  Method  of  Solution 

Coneldar  a  ray  In  direction  (9,4)  extending  froa  the  source  point  P^ 
ro  the  point  P^  on  1^^.  The  eource  region  (Eegion  I)  is  hoaogencoua  and 
isotropic;  hence,  t***  ray  ia  a  atraight  Hue  along  the  unit  vector 
First,  the  diatan.:e  r^^^  la  found  am^  the  coordinates  of  point  P^^  are 
detemlned.  Then  the  unit  vector  norml  to  the  surface  1^^  at  point  P^ 

A  A 

Is  found  (Figure  5) .  The  plane  of  vectors  r^^  and  astabllahea  the 
Incident  plane.  The  angle  between  these  two  vectors  Is  the  Incident 
sngle  o^.  Using  Snell's  law,  the  refraction  angle  Is  obtained,  which 

A 

establishes  the  direction  of  the  trenstsitted  wave,  r^2>  Seglon  11 

(dielectric) .  The  ray  In  Region  11  is  a  sti tight  line  along  the  unit 

*  *1*1' 
vector  r^.,.  Three  coordinate  syateu  >  (Uj^,v^,Mj^) ,  and 

*t  *t  * 

^*1*^1  *^12^’  with  coMon  origin  at  poiut  tre  then  established. 

They  belong  to  the  Incident  ray,  the  surface  L^,  end  the  tranaaltted 
ray,  respectively. 

The  Incident  field  (2^,^^)  Is  spilt  Into  a  noraally  polarised  field 
(ij®,5j*'),  (E  -vector  nornal  to  the  incidant  plane  at  F^)  and  a  parallsl 
polarised  field  ({^*’,3^'*).  The  transaitted  field  at  point  P^  is  obtained 
as  folloirs: 


•  r. 


1  “l 


8 


tjp  -  z 


(3.1) 


««  .  tj 


^  r 


12 


in  which  t"  Aod  nr*  trnnMinsion  confficinnts  tor  th«  no 
paraliol  polariaad  Halda.  raapactivaly^ 


I  and 


.n  2 

“rr: 


Y  cos  a 


t 


cos  o. 


nl 


"  irrj 


pi 


'nl 


pi 


Yq  cos  Oj 


cos 

_ i!!l.i 

Zq  COS  o^  ®  cos 


Z  cos  a. 


(3.2) 


Y  M  —  •  •  — i- —  Y  "  —  •  /~  •  nY  .  n  •  •  /c~  . 

^0  2q  Juq  120»  •  ^  1  "^0  •  y*o  ^ 

Hots  that  tha  aubscripc  1,  in  for  aaanpls,  signifisa  tha  fisld 
avaluatad  at  point  1. 

The  transmittad  fiald  at  P^  is  incident  on  £,  at  point  P^.  Coordioatas 
of  this  point  cau  ba  found  frost  tha  knowladga  of  the  coordinacas  of  point 

A 

P,  and  tha  transolttad  rav  direction  r,..  Tha  field  values,  in  going  fron 
X 

Pj  to  Pj.  undergo  aoaa  change  which  is  dependant  on  the  divergonca  of  tha 
ray  Thus .  wa  have 

-jkr, 


ij  •  <orj2^* 


(3.3) 


in  which  k  ■  nk^  is  tha  wave  ntsiber  in  the  dielectric  and  l)Pj^2 
divergence  factor  for  the  pencil  of  rays  travailing  fron  Pj^  to  P2  in  the 
dielectric.  It  is  given  in  (3),  as 


nr  -  n  4.  4.  H.  ^-l/2 

DF12  (1  rj^j)  (1  ^2  *12^ 


(3.4) 
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in  which  and  q2  *'*  prioclpnl  CMrvntur«i  (inwwran  of  tho  radii 
of  curvoturo)  for  tho  roy  poncil  in  Rogicn  XI.  Thoy  or#  found  fro*  tho 
curvoturo  notrlx  of  tho  tronsnlttod  roy  ot  point  P^*  aotriji  itoolf 

ia  found  fron  o  foroulo  involvint  tho  curvoturo  notrlx  of  tho  Incidont 
roy  and  that  of  tho  turfoco  point  P^.  Tho  curvoturo  notrls  of 

tho  trononittod  roy  poncil  ot  point  P^  la  oloo  uood  to  find  tho  curvnturo 
notris  of  tho  roy  ot  point  P2  incidont  upon  tho  ourfoco  £2* 

Having  tho  fiald  incidont  upon  I2  at  point  ?2*  diroction  r22« 

and  its  curvoturo  notrix,  wo  can  procood,  in  0  nonnor  ainilar  to  tho 
tronamiooion  through  to  find  tho  fit  Id  trononittod  through  £2  at  P2 
(Piguro  6).  Thuo,  0  unit  vector  Il2t  nomol  to  £2  io  obtoinod,  and 

A 

togothor  with  r^2  dofinoo  tho  incidanco  piano  ot  point  P2.  Tho  incidonco 
onglo  02  (era  *2  "  **2  *  *^12^  Again»  Sooll'a  law  ia 

invoked  to  find  tho  refraction  onglo  ^2*  angle  opocifioo  the 

roy  direction  r2^  in  Region  II.  (ou'-oid*  the  redone).  Thr^o  coerdinoto 
*  ...  "t  *t  * 

ayotono  (X2,y2ir^2^’  <» >.y2*’^;^3)  co»oo  origin  e: 

point  P2  are  then  introducaJ. 

Tho  field  ircloost  upon  I2  at  Pj  ia  roaolvod  into  porollol 

and  nomolly  polarised  fiJldr,  fron  which  tho  trononittod  fialdo  ore  found 
Ofl  followa: 


|tn  ,  ,n  |ln 

“2  2  '2 


V23 


t** 

'2  "2 


ri? 


*0'‘2 


23 


(3.5) 
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in  uhlch 


t"  - - 2_ 

'2  i  V 


n2 


'n2 


n  1 

con  a2 


'2  “  TTv 


P2 


P2 


CO*  a. 


CO*  a. 


(3.6) 


Tha  flald  at  obaarvation  point  P^  found  fron  ch*  crananlttad  fiald 


at  P2t  Buch  that, 


ij  -  (or„)  . 


■^*^0*^23 


(3.7*) 


in  which 


OP 


III 


-1/2, 


III 


23 


.  (I><,rr23)  *'‘(l  Mr^23) 


.-1/2 


(3.7b) 


and  and  av*  cha  principal  eurvaturaa  of  th*  ray  pancil  in 
Ragion  111.  Thav  *r*  obcainad  fron  ch*  curvatur*  natrix  of  Ch* 
crananiccad  ray  at  point  P^.  Thia  natrix  ia  obcainad  fron  a  fomul* 
alraady  nantionad  in  connaction  with  trananiaaion  through  1 2*  Por 
a  typical  factor  in  Bq*.  (3.6)  and  (3.7b),  Cha  following  aquara  root 
convantlon  ia  uaad: 

' 

•♦'If I,  if  f  la  raal 

f  ■  1//1  qr  •  ■ 

■♦■j|f|,  if  f  ia  inaglnary 


(3.8) 


It  should  ba  nantionad  hara  chat  wa  hava  ignored  sniltlpla  raflactions 
in  Che  dialectric  radona  throughout  our  analyaia.  Dataila  of  cha  analysis 
nay  ba  found  In  [S]. 

We  now  sunnarisa  the  final  results  obcainad  so  f.tr.  the  point 
source  at  (Figura  4)  radiates  a  spherical  wava  described  by 


iSr.e,#)  -  (Ke.f)e  ♦  q(0.»>#). 

.  1  ■■■  .  ((p  eon^eeoa  -  Q  •!!»♦)«  ♦  (F  •iaacMF  ♦  Q  co*«)r  -  F  sloea] 

0^  (3.9») 


5^  -  YqI  » 


(3.9b) 


viMra  (tfOta)  tr*  sptMrlcal  coordlQ*t««  with  oriflB  at  Pq.  Th*  patcwra 

fuDCtions  F<6,<))  *nd  Q(8,a)  tn  (3.9)  ar*  glwan.  At  polat  (Flfura  s), 

»1 

w«  daooapeaa  tha  flald  Into  two  coag^^MMota  In  tha  dlractlooa  of  (Kj^.y^),  i.a. 

i; .  .  sj  -  voi  •  M  ■  <’•“> 


At  cha  obaarratioo  |K>tnt  Fj,  wa  axpraaa  tM  flald  aa  followa 
15  -  (i5-i;  >  (i5*yl)y5  •  ^  "  *3  • 


(3.11) 


Xha  two  oofipooaoca  of  in  (3.11)  are  found  fro*  tha  ■attl*  equation 


r**' 

®3  *2 


-  (DF^j)(0Fj3)a 


-J‘S)<“12^'23> 


t;t;(aJ**J)  ntjt;(*5-y;)l  [eJ-xJ] 


r^t;t;(Jj-*5)  tjtjcyj-yj)  I  iBj-yJ 


(3.12a) 


or  aora  coapaccly, 


-  -■"'o'-'U  *  '23>,  I 

E-  •  (DF)  a  T  Ej^  . 


(3.12b) 
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tn  (3.12),  B  !•  ctM  rafraccloa  ladas  of  ch«  <li«lBcerlc,  lo  tho 
froB-opaco  wbvb  ati^or  aad  t*.  t^,  tj,  or*  th*  Boraal  oiMl  paraXlBl 
trBDMUiBoioB  cootfieloatt  at  potato  Fj.  and  F2'  taopactlvvXy ,  at  flvan 
tn  (3.2)  and  (3.6).  TYta  two  dl/argaaca  (actora  art  ttvao  Xd  (3.4)  and 
(3.7)i).  Thatr  ealeulacioaa  ccnaclcuca  cha  najor  affort  of  cha  pvaaant 
aoluctoa. 


IV.  RKrUACTIOM  AT  A  CUHVTO  DIELECTHIC  IMTEEFACB 


A .  Introduction 

Tho  rofioction  «t  a  dialactric  intarfact  ifi  of  fundaiMntal  laportanca 
in  alacCTOMgnatic  thaory.  If  tha  Intarfnca  is  uibltrarily  cur*  tha 
only  avallabla  solutiou  la  tha  orr  tltiri''ad  by  tha  gaoMtrlcal  optica  theory 
(CO).  Surh  a  aolutlon  ccnalata  of  tvo  aaln  Ingradlanta:  tha  weil-kaoim 
Fiaancl  forvulca  for  tha  tranamlaalon  and  rafiactlon  cnafflclanta  (dua  to 
A.  J.  Fraonal  In  1623)  and  a  ao-callad  “dlvarganca  factor  DF."  Surprialngly , 
tha  aclvi'  ion  of  DF  waa  darlvad  aa  f'arly  aa  1915  by  A.  Guliatrand  (6),  but 
Its  application  vaa  not  wldaly  racognlsad  In  tha  alactrotaagnatlc/optlcal 
cosnunlty  unti.!  vary  racantly.  In  19'’2,  Daacha«spa  (31,  (4)  radrrlvad 
GuUstrand'j  result  by  using  "curvatura  Hatilcas"  for  daacrlbing  curved 
aurfacas/wavafronta,  thus  raaultlng  In  greater  clarity  and  alaplar 
rnapucatloua . 

In  this  chapter,  wa  aupplaaant  Daachaaipa '  rrauics  by  giving 
explicit  forvul&a  for  calculating  varloua  curvatura  atatrlcaa  and  by 
llluatratlng  the  physical  slgniflcanca  of  DF  via  analytical  and  nuEMrical 
exaaplea.  Another  Botlvitlon  for  the  present  work  la  to  coapara  our 
aolutlon  with  tha  one  daacribad  by  SnyJar  and  Love  (7]- [9]  for  tha  sane 
problea.  It  la  shown  th«t  these  tvo  solutions  arv  not  in  agtaeaent. 

Wa  believe  that  the  Snydar-Lova  solution  is  Incorrect. 
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B,  Final  Solution  for  the  Refracted  Fields 


We  begin  with  a  statement  of  the  problem.  Two  Infinite  dielectric 
media  with  refraction  Indices  n^  and  n2  ar^t  separated  by  a  curved  Inter¬ 
face  E  (Figure  7),  which  la  described  by 

Z;  2  -  f(Xi,y)  .  (4.1) 


The  origin  of  the  (x,y,z)  coordinates  is  at  the  source  point  0  In  medium 
1.  The  source  emits  a  spherical  wave,  whose  electric  field  at  an  obser¬ 
vation  point  r  ■  (r,e,4i)  is  given  by  [for  exp(J(i)t)  time  convention! 

-jkir 


e‘<?)  .  i- 


[ep(e,4))  +  dQ(e,(^)] 


(4.2) 


where  ■  2Tf/Xj^  -  nj^u/c,  and  (r,0,4i)  are  the  spherical  coordinates 
with  origin  at  C.  The  problem  at  hand  Is  to  find  the  transmitted  field 
E  at  a  typical  point  2  in  medium  2,  and  the  reflected  field  E  at  a 
typical  point  3  in  medium  1. 

We  attack  the  problem  by  the  geometrical  optics  theory  (GO)  [3],  [4]. 
Referring  to  Figure  7,  let  us  concentrate  on  a  typical  Incident  ray  in  the 
direction  of  emanating  from  the  source  at  0.  The  "outward"  normal  to 
surface  Z  at  the  refraction  point  1  Is  N.  The  plane  defined  by  the  ray 
01  and  N  Is  the  plane  of  Incidence.  With  respect  to  this  plane,  we  resolve 
the  Incident  f  li>'ld  E^  Into  two  components :  perpend  Icular  component  E^ 


and  parallel  component  E| 

,1 


1 

II 


We  Introduce  a  scalar  u  such  that 


E2^  ,  for  perpendicular  polarization 


"II 


for  parallel  polarization 


(4.3) 


Similar  decompositions  and  notations  apply  to  E  and  E  .  Then,  the  final 


solution  derived  from  GO  has  the  following  form 


Figjra  7.  Refraction  at  a  curved  dielectric  Interface  Z 


x> 


u^{2)  -  (DF)2  T  e  ^  u^(l) 

u’^O)  -  (DF)^  Re  ^  u^(l)  . 


(4.4«) 


(4.4b) 


The  various  factors  in  (4.4)  are  explained  below:  T  and  R  are  the  well- 
known  Fresnel's  transmission  and  reflection  coefficients  (for  a  planar 
Interface) ,  given  by 


1  +  Y  * 


1  -  Y 
1  +  Y 


(4.5) 


where 


n(ccs  a^/cos  a^)  ,  for  perpendicular  polarization 


l^n  (cos  02/003  a^)  ,  for  parallel  polarization 


n  ■  (n2/nj^)  ■  relative  refraction  Index. 


The  Incident  angle  o^  and  transmitted  angle  02  related  by  the  Snell's  law 


sin  a_  ■  —  sin  e, 
2  n  1 


(4.6) 


For  n  <  1,  a  critical  Incident  angle  exists,  where 


sin  o  ■  n  ,  if  n  <  1 
c 


(4.7) 


It  >  o^,  02  defined  in  (4.6)  becomes  comply  and  the  simple  ray  picture 
shown  In  Figure  7  Is  lost.  It  Is  not  Ivmediately  clear  hov  the  present  CO 
solution  must  be  modified.  Therefore,  In  this  paper,  we  exclude  the  case 
°*1  ^  n  <  ).  The  factor  (DF)2  in  (4.4a)  is  the  so-called  "divergence 

factor"  13)  of  the  trausmltted  ray  pencil  at  point  2  In  reference  tc  point  1, 
It  Is  given  by 


1 


I 


(4.8) 


(DF)^  -  - -  —  . 

/I  +  (b7R2j, )  /I  +  (b/R22> 

Here  (82]^>822)  the  two  principal  radii  of  curvature  of  the  trannilttod 

wavefront  passing  through  point  1.  The  sign  convention  of  R2]^(<’'  ^22^ 

as  follows:  R2j^  is  positive  if  the  trananltted  rays  in  the  corresponding 

noraal  section  are  divergent^  and  R2^  is  negative  if  the  transmitted  rays 

are  convergent.  The  square  roots  in  (4.8)  take  either  positive  real  or 

negative  imaginary  value.  Thus,  (DF)2  is  positive  real  (no  focus  between 

points  1  and  2  on  the  transmitted  ray) ,  positive  .Imaginary  (one  focus 

between  1  and  2) ,  or  negative  real  (two  foci  between  1  and  2) .  The  factor 

(DF)^  in  (4.4b)  is  the  divergence  factor  of  the  reflected  ray  pencil  at 

point  3  in  reference  to  point  1.  It  is  given  by 

(DF),  «  - - ^  - -  -  - - .  (4.9) 

A  +  (c/Rj^)'  A  +  (C/R32) 

The  determination  of  the  four  principal  radii  of  curvature  (^2] 

is  the  key  to  the  present  problem.  In  Section  IV. C,  we  give  an  explicit, 

step-by-step  description  of  their  determination. 

In  summary,  for  the  refraction  problem  in  Figure  7,  the  final  solutions 
for  the  fields  of  the  transmitted  and  reflected  rays  are  given  in  (4.4). 

This  solution  is  based  on  GO.  It  is  valid  for  high  frequencies,  and  for 
all  cases,  except  when  total  reflection  occurs  (n  <  1  and 
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C.  Calculation  of  CurvturiiM  of  Refracted  Wftv>gront» 

For  an  arbitrary  lutarfaca  I  and  an  arbitrary  Incidanc  ray  01 
(Fi|ura  7),  the  calculation  of  tha  four  radii  of  curvaturac  *^2*^31 '^32^ 
is  not  a  ainpla  taak.  In  chia  aaetion,  w*  prasant  a  syutanatir  and  axplicit 
procadure  for  doing  thia  calculation. 

Coordinate  ayatawe  at  point  1.  Consider  a  ray  leaving  tha  source  at  0 
in  the  direction  (9,4),  which  intavaacta  tha  surface  I  daacrlbad  in  (A.l) 
at  point  1.  Tha  distance  a  is  detaminad  fron  tha  non-linaar  equation 

a  cos  9  «  f(x  >  a  «in  6  cos  4,  y  *  a  sin  8  nln  4)  .  (4.10) 

The  unit  vector  in  the  direction  of  the  incident  ray  is 

•  X  sin  0  cos  4  +  y  sin  8  sin  4  +  z  cos  8  .  (4.11) 

The  unit  nomal  N  of  surface  I  at  point  1  is 

N  -  ^  (-f^  i  -  fy  y  +  i)  (4.12) 

2  2  1/2 

where  4  «  +(1  +  •♦"  f  )  >  *Qd  f  *  iov  Is  the  partial  derivative 

X  y  X 

of  f(x,y)  with  respect  Co  x.  By  definioc  4  positive,  we  have  chosen  N  in 

a  A 

Che  forward  direction  with  respect  to  the  incident  tay.  Vectors  and  N 
define  Che  plane  of  incidence.  <\t  point  7,  we  introduce  four  orthononal 
base  vactors: 

A  *  ^ 

^*1*^1*®!^  for  the  incident  ray  01 
(X2»y2<^2^  for  the  tranawlttad  ray  12 

A  n  A 

(Xj,yj,tj)  for  the  reflectad  ray  13 
(u,v,H)  for  the  surface  I 
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U«  chooa* 


(4.13) 


■  «  unit  vector  notMl  to  the  plane  of  incidence 


Then  it  follow* 

u-vxN  ,  "y  for  n  »  1,2,3  .  (4.J4) 

n  n  n 

Tile  incident,  tranenitted,  and  reflected  rays  are  respectively  in  the 
direction* 

A  A  A 

<■  u  tin  -f  H  cos 

A  A  A 

*2  •  u  sin  Oj  N  coe  Oj 

A  A  A 

*2  ■  u  sin  -  S  cea 

where 

sin  Oj  •  i  ®i»  ®2  -  '  (4.15d) 

Note  that,  bccausw  of  the  particular  choice  in  (4.13),  both  and  are 
always  positive,  and  have  values  in  (0,w/2). 

Curvature  natria  of  surface  1,  At  point  1  on  surface  C,  the  followlns 
two  vectors  lie  in  the  tangent  plane  of  the  surface: 


(4.15a) 

(4.15b) 

(4.15c) 


*  (4.16a) 

^ly  -  y  +  fy  3  (4.16b) 
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wh«r«  (x,y,g)  ar*  avaluatad  at  point  1.  With  raapact  to  thu  baaa  vactora 
tha  curvaturo  matrix  of  C  la  glvan  by  [10] 


•  ,G,  - 

f,r. 

f.B.  -  a,P/ 

1  1 

1  1 

11  11 

1 

,2 

(4.17) 

A 

hh  - 

*1^ 

•l®l  ‘  ^1^1 

whara 

A  - 

**•<1  «*• 

2  2 

f  f 

!jl/2 

*  > 

I 

L  - 

1  -f  f^ 

-  f  f  ,  C- 

-  1 

f' 

1 

X 

*  1 

y 

c 

1,  ■ 

-a“^  f 

f,  -  -a"^  f 

St 

1 

1 

> 

1 

• 

1 

XX  ' 

i  xy 

•  *1 

yy 

All  (x,y,a)'a 

are  evaluatad  at  point  1.  Now  wo 

tranafar 

tha  curvature 

aatrlx  %rlth  raapact 

to  to  that  %d.th  raapact  to 

(u.v),  namely. 

o 

t-i 

■ 

A-^ 

Qj.  A 

(4.18) 

whara 

’^Ix 

•u 

r,  “V 
lx 

A  - 

• 

♦ 

• 

•u 

r,  ‘v 

ly  J 

It  awy  ba  ahown  (10]  that  a  principal  radlua  calculated  from  (4.17)  or  (4.18) 
haa  a  poaltlva  (negatlva)  algn  If  tha  normal  aactlon  of  tha  aurfaca  banda 
away  from  (toward)  tha  normal  N.  For  axampla.  If  I  la  a  aphara  with  radlua 
p  and  tha  normal  N  polnta  away  from  tha  aphara  cantar,  wa  hawa 


(4.19) 
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U*  noc«  that  th«  praaaot  algn  cottv«nwlon  for  tho  ourfaco  curvatura  la  tha 
aoM  aa  chat  uaad  la  (4],  but  oppoaita  to  chat  in  (3),  (101. 

Curvatura  matrlcaa  of  <wiaafronta.  Tha  Incldanc  uo'/afroat  paaolug 
through  point  1  la  apharlcal  with  a  radiua  a.  Thua,  Ita  curvatura  natrix 
with  raapact  to  baaa  vactora  or  any  othar  orthonoraal  baaa 

vactora,  is 


(4.20) 


Tha  curvatura  natricaa  of  tha  trananittad  and  raflactad  wavafronta  paaalng 
through  1  ara  axpraaaad  with  raapact  to  baaa  vactora  (x.tty2)t  and  (x^ty^), 
raapactivalv .  Thay  ara  danotad  by  and  Q^.  Tha  solution  of  Q2  ia  found 
fron  tho  following  natrix  aquation  (4): 


n  Bj  •  ^  (n  coa  *2  ”  c®*  Qj;  (4.21) 

whara 


%  A  a  A 

X  U  Xj^*v 

com  0 

*1  •  * 

m 

L>1'U  y^.v 

0  1 

“a  A  A  a' 

cos  82  0 

■ 

m 

£ 

A  A  A  A 

yj’u  yj'v 

.  0  K 

The  solution  of  ia  found  fra«  tha  following  natrix  aquation 
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(4.22) 


Ij  Qj  Bj  •  bJ  B^  -  2(coa  a^)  Qj. 

whara 

•coa  0 

0  1 

Principal  radii  of  curvatuta  of  rafractad  wavafronta.  Onca  aacricat 
Q2  and  ara  dacarainad  front  (4.21)  and  (4.22),  chay  nay  ba  diafonaliaad 
in  a  acandard  nannar  Co  find  chair  aiganvaccora  (principal  diracciona  of 
Cha  wavafronc)  and  chair  aiganvaluaa  (principal  curvacuraa)  [10).  In 
particular,  cha  principal  radii  of  Cha  tranaaiccad  twvafronc  (^21*^22^ 
cha  rooca  of  cha  following  quadratic  aquation 

^  ^  (craca  Q^)  dac.  Q2  •  0  .  (4.23) 

R 

If  1)2  (4.23)  la  raplacad  by  Q^,  Cha  cvo  rooca  ara  tha  radii 

of  cha  rafleccad  vavafronc. 


*3  • 


Xj‘u 


yj*u  yj*w 
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D.  Sp<ci>l  SphTical  Infrf«c« 


To  llluotrat*  th«  rMult*  obe.olnod  in  th«  pravious  two  ••ctiona, 
l«t  us  concsatrata  on  s  spaclsl  csss  in  which  ths  intsrfscs  £  is  sphsrirsl 
with  rsdius  |p|t  ss  shown  in  Figurs  6.  Following  our  sign  convention,  ths 
rsdius  of  curvsturs  of  I  in 

,  if  £  is  conesvs  when  viswsd  fro*  tha  source  (Figure  6s) 


if  £  is  convex  (Figure  8b) 


(4.24) 


Without  loss  of  geners7.ity,  we  sssune  thst  the  incident  rsy  fron  the  source 
at  point  0  is  In  the  direction  (9,d  *  0).  The  plane  of  incidence  is  than 
the  X  >  s  plana.  Making  use  of  the  foraulss  in  Section  IV. C,  we  find  thst 
the  principal  radii  of  the  ttansaitted  and  reflected  wavefronts  passing 
through  point  I  are 


■  (n  cos 


“2>  [l 


2  1  , 

cos  a,  +  -  (n  cos  a_  -  cos  a 
10  Z 


.>] 


-I 
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1  a-  I  X 

—  +  —  (n  cos  a-  -  cos  o, ) 
na  no  2  1 


-1 


-  ^ 

[a  p  cos  a^J 

.  [i . 

L*  °  J 


(4.25e) 

(4.25b) 

(4.26a) 

(4.26b) 


It  can  be  shown  that  R2^  and  R^^  tadii  of  curvature  of  the  nomal 

sections  in  tha  x  -  s  plane  (plane  of  incidence) ,  whereas  R22  Md  R^2 
arc  those  in  the  orthogonal  directions.  Since  in  general  R2^  f*  R22 
*^32*  the  refracted  and  reflected  pencils  are  astigaatic. 
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(A)  Wor—1  Incldtnca.  For  •  0 


(^.2))  and  (4,26)  b«co«« 

(4.27) 

(4.28) 


Thus,  for  norMl  incldmc'i,  both  rofractod  poncili  h4va  •pharirol  vavofronts 
(no  longar  •■tigaaclc).  Th«  rtilaeion  in  (4.27)  aay  ba  rearranged  co  taad 


(4.29) 


which  Is  tha  vail>kr.own  lana  equation  in  optics.  t^Sse  for  example 
Eq.  (40-14),  p.  678  of  [11].)  Note  the  correapondinf  notatiotia  uaad  in  [11] 
and  ^ere:  n  -*  1,  n'  -*■  n,  a  -►  «,  a*  -►  (“*^21^’  ^  *  (“O).)  The  divergent 

iacid  'nt  pencil  from  a  point  sour>'a  ia  converted  li.Cu  e  convergent  trenaaitted 
penci.1  in  Medium  2  %rt»en  <  0.  This  occurs  when 


nj 


>  0 


(4.30) 


If  0  >  0  (concave  dielectric  intorfece  shown  in  Figure  8a),  this  is  poauibla 
if  n  ■  (nj/Uj^)  <  i.  If  0  <  0  (convex  dielectric  interface  ahovn  in  Figure  8b), 
this  is  possible  if  n  >  1. 

(B)  CoMparis^n  with  Synder  and  Love’s  reeult.  In  a  recant  article  [7], 
Synder  and  Love  consider  the  problem  sketched  in  Figure  7  for  an  incident 
plane  wave  (source  distance  a  in  Figure  7).  Their  final  result  is  in 
disagreement  with  ours.  To  show  this  diaagrcrment,  let  ua  concentxate  on 
a  simple  case  (Figure  9):  a  concave,  spherical,  dielectric  interface  ia 
illuminated  by  an  incident  plane  wave  %d'.,..ch  is  given  by 
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1  - 
E  (x,y,t)  -  I  • 


(4.31) 


TViC  problaa  is  co  fin;)  the  high-fraquancy  iioluClon  of  cha  refractad  lialda 
along  tha  yaxia.  lUaad  on  CO,  our  solution  la  givan  In  (4,4),  (4.5),  (4.8), 
(4.9),  (4.27),  and  (4.28}  with  »  *  Using  tha  coordinata  systaa  in  Figura  9, 
tha  transaitted  and  raflectad  fialda  ara 


r(0.y,0)  •  ^ 


1  +  n 


e'<0.,.0)  .  «I_1  1  i2|y|/„-J  Hr  •  ‘ 


V  >  0 


y  <  0 


(4.32a) 


(4.32h) 


Tha  factors  in  [  ]  in  (4.32)  a^e  diverganca  factoiu.  The  intensity  or  power 
denaiC'/  of  Cha  incident  field  on  tha  central  ray  (along  x  ••  t  ■  0)  la  given  by 


■  he{y(E^  *  H^*) }  ■  (nj^/120r)  wates/a^ 


(4.33) 


which  is  independent  of  y.  This  is  due  to  rhe  fact  chat  the  incident  field 
ia  a  plane  wave  and  ell  incident  rays  are  parallel.  Tha  intenaity  of  the 
refracted  field  on  the  centriil  ray  does  vary  with  y,  namely. 


iMv') 


_2_ 

1  +  nj 


1  t  |2^1  ^^1 

{  n  J  IpJ 


,  y  >  0 


(4. J4a) 


[1  -  nl 

2 

1 

[i  +  nj 

.1  -  (2|y|/p). 

y  <  0 


(4.34b) 


At  the  focal  point  of  the  reflected  pencil  v  •  -p/2  in  medium  1,  the  intensity 
in  (4.i4b)  predicted  by  Che  present  GO  becomes  infinite  as  expected.  For 
the  same  problem  sketched  in  Figure  9,  Snyder  and  Love's  solution  is  given 
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in  Fqs.  f29)  through  ot  [7].  For  the  case  of  noraal  Incidence 

(a^  ■  0)  and  central  ray  (x  ■  z  -  0),  Snyder  and  Love's  solution  reads 


^SL(O.y.O)  -  z  .  p  »  y  >  0 

TCiT  1  -  n  ■*  1  .  /s 

?s)/o.y.O)  -  z 


(4.35a) 

(4.35b) 


which  should  be  compared  with  our  solution  In  (4.32).  We  note  that 
(i)  dive/ge-xe  factors  (DF)2  and  (DF)^  are  missing  in  (4.35),  and  (11)  the 
propagation  phase  factor  exp<^-jk2y)  is  missing  In  (4.35a).  Thus,  we  believe 
that  (4.35)  is  Incorrect.  Furthermore,  for  each  Incident  ray  (fixed  a^). 
Snyder  and  Love  define  a  "power  transmission  coefficient  by 

(Eq.  (35a)  of  [7]) 

f  ^ 


SL 


(4.36) 


As  may  be  seen  from  (4.34h),  the  intensity  I  is,  in  general,  a  function 
« 

of  position  (x,y,z)(  because  of  the  divergence/convergence  of  the  reflected 
rav  pencil.  Then,  when  calculated  correctly  is  also  a  function  of 
position,  and  does  not  have  the  usual  significance  assoc:'.ated  with  the 
"power  tr '.nsmieeion  coefficient." 
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K.  tlumerict.1  Reaulta  «ad  Dtscugsiot;. 

For  Lhe  refraction  probleji  sketched  in  Figure  the  flnftl  solution 
for  the  transaltted  field  into  aedlun  2  Is  glvei^  lr<  (A. 4s),  when 
the  Incident  field  is  gi''en  in  (A.t).  In  this  section,  we  present  aome 
numerical  results  for  for  various  intixrfscas  and  source  locations. 

We  consider  three  types  of  interfaces:  the  sphevlcal  Incorfacc 
described  by 

(r/X^)  -  i  -  [1  -  (x^  +  ,  ,4.37) 

the  paraboloidal  Interface  described  by 

iz/\^)  •  (x^  +  y^)/2xj  ,  (4.38) 

end  Che  hyperboloidal  interface  described  by 

(z/Xi)  -  i(l  +  2(x^  y^)/xjj^^^  -  I  (4.39) 

where  X^  is  Che  wavelength  in  medium  1  in  wiiich  the  source  is  located , 

For  easy  r.oavarlson,  we  have  chosen  the  above  interfaces  such  that  they 
all  have  tne  -^aae  curvature  in  the  axial  direction  (Figure  1C).  There  are 
six  source  locations,  indicated  by  numerals  inside  a  small  circle.  When 
the  source  is  at  location  1,  for  example,  medium  1  is  on  the  right  and 
medlua  2  on  che  left,  and  the  Interface  is  concave.  The  source  is  assumed 
to  be  y-polarlzed.  We  calculate  the  transmitted  field  in  the  E-plane 

A  A 

(plane  normal  to  x)  and  H-plane  (plane  normal  to  y) ,  In  these  two  planes, 
Che  Incident  field  is  assumed  to  be 


2.0 


KYPERBOLOIO  PARABOLOID 


Figure  10.  Three  dielectric  loCerfeces.  At  x  •  0,  ell  three  interfeces 
have  the  sene  rediue  of  curvature  of  1 


E-pl«ne 


(4.40) 


i^(v) 


A 

0  1 

y  1 


H-plane 


Thus,  la  Che  E-plane,  Che  f-vecCor  Is  parallel  Co  Che  plane  of  incidence; 
whereas  in  Che  H-plane,  Che  l-vecCor  is  perpendicular.  The  observacion 
poinc  2  is  in  nediua  2  (Figure  7)  wich  disCunce  b  -*-  «•  (far  xone) . 

We  calculaCe  Che  normalised  far  field  defined  by 


1^(2) 

E- field  when  n^  ^  n2 

t^(2) 

E-fleld  when  02*02 

SubsCiCuCa  (4.4a)  and  (4.40)  inco  (4.41).  Under  Che  condicion  b  ->■ 


(4.41) 


we 


have 


EN 


b  • 


(4.42) 


Here  a  is  the  disCance  between  Che  source  and  Che  inCerface  along  Che 
incidenC  ray,  and  T  is  Che  Fresnel's  Cransmisuion  coefficienC  given  in  (4.5). 
The  facCor  *^^21^22  radius  of  Che  Gaussian  curvature.  In  presenting 

Che  numerical  results,  we  plot  EN  as  a  function  of  0,  where  0  is  Che  polar 
angle  of  observacion  point  2  measured  from  a  line  parallel  Co  Che  z-axis 
and  passing  through  the  source  point.  The  relative  index  n  ~  **2^^1 
is  always  set  at  2  (cransmlasion  into  a  denser  me  Hum). 

Concave  spherical  Interface.  Figures  11  and  12  show  the  E-  and  H-plane 
.ar-fleld  pattern  EN  as  a  function  of  0.  Note  Chat  Che  field  strengths 
Increase  as  the  source  moves  closer  to  the  Interface  (smaller  a) .  This 
is  mainly  due  to  the  face  that  EN  is  inversely  ; roportional  Co  a,  according 
Co  (4.41).  The  Gaussian  curvature  <^21*22  <1**^^*****  wich  a„  but  not  enough 
to  offset  Che  factor  (l/a)  in  (4.41).  For  source  3,  which  is  at  the  center 
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POLAR  ANGLE.  0(DE6) 


of  th«  ophorical  intarfaco,  all  of  eha  incldant  raya  ara  nonal  to  tha 
latarfaca.  It  can  oa  ahown  chat  II21  *  *^2  *  **  ^  calculated 

fron  (4.41)  la  equal  to  which  la  0.667  for  tha  praaant  caaa  of  n  ■  2. 

Of  particular  Intaraat  la  tha  H-plana  pattam  of  aourca  1  ahoim  In  Fltura  12. 
Note  tha  aarkad  aaynMCry  In  tha  far>flald  pattern  which  la  due  to  the 
aayMatry  of  tha  aurfaca  with  raapact  to  aourca  1.  Figure  13  ahowa  tha 
variation  of  the  axial  far  field  whan  tha  aourca  la  novad  along  and 
parallel  to  tha  i-axla.  It  ahowa  clearly  tha  Incraaaa  of  the  field  aa 
tha  aourca  novas  closer  to  the  Interface. 

Concave  paraboloidal  Interface  (Flguraa  14  and  15) .  Tha  far-flald 
pattama  for  the  concave  paraboloid  ara  quite  slnllar  to  those  for  tha 
spherical  case.  However,  tha  pattern  variations  ara  nora  pronounced. 

Concave  hyperbololdal  Interface  (Figures  16  and  17).  Note  that  the 
far-flald  pattam  due  to  source  4  has  a  dip  Instead  of  a  peak  In  the 
axial  direction.  This  la  In  contrast  to  the  situations  In  Figures  11,  12, 

14,  and  15.  There  Is  another  fact  worth  awntlonlng.  Because  of  the  choice 
of  tha  sane  axial  curvature  for  the  above  three  interfaces,  the  axial  field 
Is  the  sane  for  al?.  Interfaces  when  the  source  is  at  2,  3,  or  4.  However, 
for  source  1,  which  la  displaced  fron  the  synnetry  axis,  the  nomallsed 
axial  field  EII(9  >■  0)  Increases  fron  0.826  for  the  spherical  surface  to 
0.954  for  the  hyperboloid. 

Convex  Interfaces.  Tha  B-  and  H-plana  far-flald  patterns  for  a 
convex  sphere,  paraboloid  and  hyperboloid  are  shown  In  Figures  18  and  19 
for  source  locations  5  and  6.  The  source  locations  5  and  6  were  chosen 
based  on  (4.30).  Source  5  produces  a  divergent  axial  pencil  In  nedlun  2, 
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-field  puttern  through  •  concave  paraboloid 


Figure  15.  S«W!  an  Figure  14,  except  for  B-plane  pettem. 


bolold 


POLAR  ANGLE,  a (0£G? 

Pigure  1/.  Smk  ait  Pigrre  ifc,  cuccpt  for  H-(  ine  pattcra. 
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Figure  18.  E-plaiic  far-field  pattero  tiirough  convex  interfaces. 


Figure  19.  Saive  nn  Figure  12,  exccpc  for  M~plMoe  pactern 


whereas  source  6  produces  a  convergsnt  axial  pencil;  the  behavior  in  tli-i 
non*>axlal  direction  la  governed  by  the  type  of  the  Interface.  Thus, 
tAay  be  seen  from  figures  18  and  19,  the  far  field  in  the  axial  direction 
through  the  spherical  interface  has  a  peak,  for  source  6  and  a  dip  for 
source  5.  This  is  also  the  case  for  the  paraboloid.  However,  this 
behavior  Is  not  observed  in  the  hyperboloidal  pattern. 

For  all  the  convex  interfaces,  the  variation  of  EM  as  a  function  of 
6  In  (4.41)  is  predominantly  determined  by  the  radius  of  the  Gauaslsn 
curvature,  *^21^22*  ^  leaser  extent  by  T  or  a. 

Ray  Plctt.ro.  The  H-plane  pattern  due  to  source  6  for  a  convex 
sphere  is  given  in  Figure  19.  The  corre.spcudlng  ray  picture  is  sViown 
in  Figure  20.  We  launch  6  rays  at  4*  apart  in  the  upper  half  x-z  p3ane 
(x  >  0).  The  transmitted  rays  are  first  convergent,  and  after  crossing 
the  caustic  surface,  become  divergent.  The  Incident  rays  in  the  upper 
half  x-z  plane  within  a  20**  angle  give  rise  to  transmitted  rays  in  the 
lower  x-z  plane  (x  <  0)  within  a  13.5“  angle.  There  .ite  two  caustic 
surfaces  associated  with  the  transmitted  rays.  The  intersections  of 
the  caustic  surfaces  and  the  x^z  plane  are  indicated  by  crosses  and  dots. 
Similar  ray  pictures  can  be  drawn  for  the  other  cases  also. 
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F.  Conclualon 


lor  th*  refrmctlor.  ^roblom  sketched  in  Figure  1 ,  the  final  geooietrlcal- 
opClc<i  solutions  for  the  transmitted  field  and  the  reflected  field  are 
given  in  (A. 4).  Tiiey  are  applicable  under  rather  general  conditiona, 
namely,  tne  dielectric  interfere  described  in  (4.1)  is  arbitrary,  and  the 
incident  fisld  in  (4.2)  from  a  point  source  is  arbitrary.  A  major  step 
ir  calculating  these  so  .utio.ts  is  the  evaluation  of  the  divergence  factors 
in  (4,8)  aii'^  (h.9),  vttlch  involves  the  matrix  operation  des  cribed  by  (4.21) 
and  (4.22).  (Strictly  speaking,  the  present  solution  la  valid  in  the  high- 
frequency  limit  u)  -*■  <»;  heuevssr,  practical  experience  has  shown  that 
soLutlons  of  the  present  type  aru  reasoi-iably  accurate  as  long  as  the 
radii  of  curvuti're  of  the  dielectric  Interface  are  in  the  order  of  a 
wavelength  or  laorc. 
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V.  WAVE  TRAMSMISSIOM  IWROUGH  A  SPHEBiCAL  DIELECTRIC  SHELL 


A.  Introduction 

One  of  the  fundaaentel  problems  In  electromsgnctic  theory  Is  the 
trsnsaission  of  s  spherlcsl  wave  through  a  dielectric  shell.  This  problem 
has  numerous  applications  In  antenna  radomes,  electromagnetic  shielding,  and 
scattering.  It  appears  that  solutions  to  this  problem  are  available  only 
for  the  special  case  where  the  shell  Is  an  infinite  dielectric  half  space. 
That  case  was  first  studied  by  SoMerfeld  In  1909,  whereas  later  research 
was  summarized  in  a  book  by  Brekhovsklkh  (Chapter  IV  of  [12]).  In  this 
chapter,  we  consider  u  more  general  case,  namely,  ths  shell  has  two 
spherical  boundary  surfaces.  Unllks  the  Sommerf eld's  problem,  our  case 
does  Include  the  effects  of  the  shell's  ;-.urvature  and  thickness.  Therefore, 
Its  solution  should  be  of  more  practical  Interest. 

To  solve  our  problem  rigorously,  the  spherical  wave  expansion  may  be 
used.  However,  due  to  the  fact  that  the  source  location  and  the  two 
dielectric  surface  centers  do  not  coincide,  the  translatloual  addition 
theorem  for  vector  spherical  wave  function  [13]  must  be  used.  (Our  problot 
Is  roughly  comparable  to  scattering  by  three  dielectric  spiieres.)  This 
theorem  leads  to  a  complex  series,  tihlch  mkes  it  very  difficult  to  generate 
numerical  results.  In  this  paper,  we  use  the  geometrical  optics  theory  (Gu) 
[4]t  [14]  to  calculate  the  transmitted  field  In  the  problem  sketched  In 
Figure  21.  Such  a  solution,  though  only  approximately  valid  for  high 
frequencies.  Is  given  in  a  simple  closed  form.  Thus,  It  allows  us 
to  study  the  "cause  and  effect"  of  the  various  parameters  lu  a  convenient 
marjier  and  gain  physical  Insight. 
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B.  Solution  for  Conontric  SphTical  Sh«ll 


Lot  ut  consldor  tho  tmv*  tranaalMlon  problaa  In  Flgurt  21.  Th« 
boundnry  aurfncM  of  th«  dlnlnctric  nholl  nr*  two  concantrlc  •phono 
with  radii  and  R2  (both  poaltlva),  and  with  a  coaawn  cantor  at  Q. 

Tha  point  aourca  at  0  awlta  a  apharlcal  wava  daacrlbad  by  (for  axp  Jut 
tlaa  convantlon) 

l^(r)  -  ----  lP(e,#)e  ♦  Q(e»d)dl  .  (5.1) 

Hara,  (r,6,d)  ara  apharlcal  coordlnataa  of  r  with  origin  at  0.  Fdnctlona 
P  and  Q  doacriba  tha  radiation  pattarn  of  tha  aourca.  Tha  wava  nuabar 
Icq  •  -  u(wq£q)^^*  la  that  of  fraa  apaea.  By  ualng  GO  (4],  [14 J, 

tha  tranaalttad  flald  at  an  obaarvatlon  point  3  la  to  ba  datacalnad. 

Without  loaa  of  ganarallty,  wa  aasuaa  that  point  3  .  In  tha  (Xyi)-plana 
with  rectangular  coordinates  (x^tO^z^). 

A.  Ray  Tracing.  In  accordance  with  Snail '•  law,  we  trace  a  ray 
iron  aourca  point  0  to  obaarvatlon  point  3,  via  refraction  points  1  and 
2  (Figure  21).  Clearly,  all  four  polnta,  0  to  3,  11a  In  tha  aoM 
(x,z)-plana.  For  a  given  launching  angle  9  of  tha  ray  and  tha  dlatanca 
c,  the  other  gaooatrlcal  parawatara  can  ba  datamlnad  froa  tha  following 
rclatlona: 

aln  ■  (d/R2)aln  9  ,  a  •  R2[aln(9  -  a^)]/aln  9  (5.2) 

•In  Oj  •  n”^  aln  aj  ,  aln  o^  ■  (R^/Rj)  »ln  o^ 

^  “  02)l/oin  ,  sin  ■  n  aln 

X3  -  a  sin  9  +  b  •ln(9  -  aj  aj)  +  c  aln  (oj  -  aj  -  aj  +  6) 


49 


•  cot  0  ^  b  coa  (•  ~ 


*■  c  coa 


“2  -  - 


ttj  ®) 


Thua,  for  a  tlvan  (b.c),  wa  can  dacarvlaa  tha  poalclon  of  point  3 
atraifht forwardly.  On  tha  othar  hand,  whan  point  3  la  axpllclt 

forwulaa  do  not  axlat  for  datarmlolng  0.  Ona  haa  to  find  9  by  trial* 
and-arror. 

B.  Flald  on  tha  Bar.  Tha  praaant  vactor  flald  problM  can  ba 
dacoapoaad  Into  two  acalar  onaa:  ona  with  tha  alactrlc  field  vactor 
parpandicular  to  tha  plana  of  Ineldanca  (f  ■  y  E^)  and  tha  othar  with 
tha  alactrlc  flald  vactor  parallel  to  tha  plana  of  Ineldanca  (1||)  • 

Tha  final  aolutlon  for  tha  trananlttad  flald  at  point  3  darlvad  by 
GO  la  glvan  by 


^(3) 


e'  (3) 

eJ  (1) 

-  (DF) 

bS  (3) 

t„  tl,  (1) 

JkQ(nb+c) 


(5.3) 


Tha  various  factors  in  (S.3)  ara  axplalnad  balow:  and  T^  ara  tha 

products  of  cha  trananisslon  coafflclants  at  points  1  and  2»  glvan  by 


r 

-1 

r  ti 

cos  a. 

1  *•> 

.  .  A  - 

A  *  j 

n  4 

COS  Cla 

aS  X 

coo  a. 

X 

(5.4a) 


T„  -  4 


ti 

-1 

r  tn 

1  ®1 

cos  0, 

■1  _  * 

X  m 

n  1 

cos 

1  ♦  n 

coo 

-1 


r5.4b) 


Tha  two  conponants  of  tha  Incident  flald  ara  calculated  fro*  (S.l): 


v(i)  - 

o 

■ 

m 

O 

cy 

-Jkn* 

E„^(l)  - 

P(e.p  -  0)  . 

(5.5a) 

(5.5b) 
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Th*  dcrlvaclon  of  th«  dlvortaneo  factor  DF  In  (5.3)  raqulraa  som  affort. 


With  cha  datailo  gtvan  in  Appandix  A,  tha  final  raault  for  DF  raada 

At  point  1  or  2,  It  can  ba  ahown  chat  cha  principal  dlractlons  of  cha 
cranaalctad  wavafronc  art  praciaaly  tha  two  dlractlons  parallel  and 
parpandleular  to  tha  plana  of  Incldaaca.  Bara  ^ 

principal  curvaturaa  of  tha  transalttad  wavafronc  at  point  1,  and 
^*^21* *^22^  ara  thoaa  at  point  2.  Thay  ara  calculatad  froa  tha  ralaClor.s 

<11  ■  (n  coa^  ®1  ■*“  “l  “  “•  *1^]  (5.7a) 


'12  *  k  *  r<“*  “i  -  ^  “*  “i> 


(5.7b) 


(coa^  *2^”^  *11^'*^  “  ®2  alj  -  n  coa  a^)] 


-lb-1 _ 2  1^1 


(5.7c) 


/V  •  “Iv”!  .  1  /  t  1.. 

<22  •  "vb  y  <12)  >  Oj  -  n  cos  Oj)  . 


(5.7d) 


Tha  sign  convention  of  x  la  as  follows.  If  x  la  poaltlva  (nagativa) , 

tha  normal  aactlon  of  tha  wavafront  is  dlvargant  (convargant) .  For 

axaapla.  If  tha  transalccad  wavafront  at  point  1  la  cha  aaaa  as  cha 

incldant  apbarlcal  wavafront,  tfa  hava  Kh  •  <12  •  -Hli^.  For  a  typical 

factor  In  (5.6),  Cha  square  root  convention  is  . 

♦!f  I ,  If  f  la  real 

f  -  (1  +  Kb)"*^^  -  ■ 

,<>j|f|.  If  f  is  laaglnary 

Whan  f  la  imaginary,  It  aaans  that  tha  ray  has  crossed  a  focus  ol  tha 


(5.8) 


ray  pencil.  Tha  (-fj)  accounts  for  tha  well-known  (ir/2)  phase  retardation. 
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Several  generel  co— inta  about  the  aolutlon  in  (S.3)  ere  in  order: 

(1)  For  the  case  whore  e  total  reflection  occura  (o^  or  becoaea  coeplax) , 
the  field  in  the  tranaaitted  region  la  not  an  optical  field,  and  the 
preaant  ray  aolutlon  (5.3)  la  no  longer  valid.  (11)  Except  for  apecial 
caaaa,  a.g.,  noraal  Incidence  •  0,  the  two  curvaturaa  (''21*'^22^ 

Che  trananltted  wavefront  aaerglng  froa  the  dielectric  ahell  are  not 
equal.  Thua,  the  tranaaitted  pencil  la  generally  aatlgaatic.  (Ill)  It 
la  poaalble  that  <2^  and/or  <^2  n*t*tlva.  Then  the  divergent  pencil 
froa  the  aource  la  tranaforaed  Into  a  convergent  (focuaing)  pencil  after 
propagating  through  the  dielectric  ahell.  (Iv)  The  aolutlon  In  (5.3) 
reaalna  valid  for  wore  geoaeCrlea  Chan  the  one  ahown  in  Figure  21.  Thla 
la  dlacusaed  further  In  the  next  section. 

C.  Generalliatlon  of  Final  Solution 

For  tranaalaalon  through  the  concentric  spherical  ahell  In  Figure  21, 

Che  final  solution  consists  of  two  parts: 

Parc  A:  Ray  tracing  foraulas  in  (5.2) 

Part  B:  Field  solution  in  (5.3)  through  (5.5) 

It  can  be  shown  Chat  Part  B  is  valid  under  a  aore  general  condition 
(Figure  22),  nanely, 

(1)  The  centers  of  the  spherical  shells  and  Q2  need  not 

coincide,  as  long  as  the  four  points  (0,1, 2,3)  are  coplanar. 

(11)  The  surfaces  of  the  shell  can  be  either  concave  or  convex. 

Looking  froa  the  aource  side,  (or  R2)  la  positive  If  the 
surface  Is  concave,  and  R^  la  negative  If  the  surface  Is 
convex. 


52 


(a)  R,  >0 ,  R2>0  (b)  R,  >0 ,  R2<0 


x-2  PLANE 


Cc)  R,  <0 .  Rt>0  (d)  R,  <0 .  Rt<0 


Figure  22.  General  conf Iguratione  where  the  field  solution  (2.3)  can  be 
used.  Points  (0.1. 2. 3)  auat  be  coplarar. 


Condition  (i)  it  ntcttstry  in  ordti  vor  th*  tctlar  trtamittion  cotfficltntt 


tnd  'V||  in  (5.4)  to  bt  valid.  Mhan  tha  four  points  (0,1, 2, 3)  art  not 
coplanar,  tha  parpandicular  and  parrllal  polariiationa  are  r.r  lonyar 
uncouplad.  Than  tha  scalar  or  Ty  in  (5.3)  aust  ba  ruplacsd  by  a 
(2  X  2)  BMcrix.  Alno,  (5.7)  baeoaas  aora  cosiplicatad. 

For  tha  ganaral  configurations  in  Plgura  22,  the  rsy  tracing 
foraulaa  in  (5.2)  ara  not  valid.  Howavar,  by  following  Snail’s  law, 
tha  ray  tracing,  avan  in  tha  aoat  ganaral  situation,  is  conceptually  siapla. 
Thus,  instead  of  working  out  a  set  of  ganaral  fonsulas,  wa  laava  it  to  the 
individual  problaaa. 

D.  Axial  Incidence  on  Syaatrical  Shells 

To  study  tha  features  of  the  present  ray  solution,  let  us  concentrate 
on  a  special  case,  where  the  four  points  (Q^,Q2,0,3)  ara  on  a  straight 
line  (Figure  23).  Then  the  four  curvatures  in  (5.7)  reduce  to 


PC  .  PC  -  ^  -v 
11  *12  no  n  Rj 


*22  “  *21 


nRj^  an(n  -  1) 
bRj^  sb(n  -  1)  na 


1  -  n 


(5.9a) 


- -  .  (5.9b) 

b  *  «!}  "2 


An  irtsrastlng  question  is  when  does  or  K22  negstive  (sManing 

a  convergent  pencil)?  This  is  anavarad  below: 

(1)  Negative  tranaaitted  pencil  inside  the  dielectric 

shell  is  a  convergent  one  when  ^  0  or 


r-^  R,  >  0 
1  -  n  1 


(5.10) 


If  Rj^  >  0  (concave  dielactric  interface  shown  in  Figure  24a),  this  Is  possible 
if  n  <  1.  If  R^  <  0  (convex  dielactric  interface  shown  in  Figure  24b), 


chit  it  pottiblt  only  If  n  >  1.  Tht  ditctnct  btcwten  point  I  tnd  focal 
point  it 

(ii)  Netativa  iC22*  Tha  cranaaitcad  fiald  in  cha  fraa-apaca  ragion 
ouctlda  tha  dlalactric  thall  it  convargant  whan  *22  <  0*  L*i  u*  concancraca 
on  a  apacial  caaa  in  which  cha  chicknaaa  of  cha  dialactric  thall  it  nail 
to  Chat 

b  «  kTJ  .  (5.11) 


Than  <22  (^•9b)  bacoaaa  approximaCaly 


'C22  -  -  (n  -  1) 


^  1_ 
R,-  Rj 


(5.12) 


which  it  Cha  wall»known  aquation  for  cha  chin  lana.  (Saa  for  axaapla 

Eq.  (41-1),  p.  665  of  [11].  Nota  cha  corratpondlng  nocacionc  utad  in  [11] 

« 

and  hart;  t  ■  -a,  t'  •  "<22*  ^  "  '*^1*  ^  *  "^2*^  condition 

for  a  nagatlva  undar  Cha  approxiaaCion  in  (5.11)  Is 

R1R2 


(n  -  l)(Rj^  -  R2) 


>  0 


(5.13) 


(ill)  Far  Fiald;  If  Cha  obtarvaclon  point  3  it  in  tha  far  lona 


(<22*^  >>  I),  than  DF  in  (5.6)  bcconas 


« 

U 

1  ^  (n  —  1  ^ 

*^2  '  *4  I 

i-i 

1  ’ 

a  +  —  b  +  c 

'■  u 

1 1 

•^1*^2 

b  n 

a 

DF 

Tha  firtc  factor  in  {  }  in  (5.14)  it  cha  divarganca  factor  of  a  planar 
dialactric  tlab  (R^,R2  **-  **) .  Thua,  cha  ratio  of  tha  alactric  fiald  at 
a  far-fiald  point  3  for  a  tpharical  ahall  and  that  for  a  dialactric 
slob  it 


-1 

K5.14) 
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_ |f^(3)  for  spherical  dhellj _ 

|e^(3)  for  a  slab  of  same  thlckaess 


1  +  (n  -  1; 


ab 

[■^2  -  ^ 

R-R- 

b 

_  1  +  1 
n 


fl  ^1 

1  -y ' 

1  ' 

a  J 

^  J 

-1 


(5.15) 


As  a  aiiuerlcal  example,  conslUer  the  case  in  which  the  luoer  and  outer 
dielectric  surfaces  are  coii.^ave  and  concentric  (Qj^  ■  Qj)  with  (Rj^/b)  *  2, 

We  plot  n  as  a  function  of  (a/b)  for  n  ■  0.5  and  a  »  3  in  Figure  25.  We 
note  that  n  can  be  substantially  different  from  unity.  When  a  >  R^,  we 
have  n  ■  1.  Thus,  In  this  interesting  special  case,  the  axial  far  field 
through  a  concentric  spherical  shell  and  that  through  a  planar  slab 
become  the  sarae.  Another  Interesting  special  case  cccurs  when  DF  ». 

It  means  that  the  paraxial  rays  emerging  from  the  dielectric  shell 
(Figure  23)  are  parallel  to  the  axis  so  that  they  focus  at  the  far-fleld 
point  at  Infinity.  From  (5.14),  It  is  shown  that  DF  ^  if  0  or 

.  R  -  R_  +  b(n  -  l)/n 

I  ■  <“  -  »  -  l.(n  -  1)/,1  • 

Under  the  thln-lens  approximation  b  I!  0,  (5.16)  is  reduced  to  the  well- 
known  lenemaker's  equation  (see  F.q.  (41-2),  p.  685  of  [11]).  In  the 
antenna  radome  appllcatloa,  (5.16)  Is  useful  In  Che  determination  of  the 
enhancement  of  the  antenna  main  beam. 

(iv)  Multiple  Refraction;  For  a  given  source  point  0  and  observation 
point  3  in  Figure  21  or  22,  we  can  trace  two  types  of  geometrical  optics  rays. 
The  first  type  Is  the  direct  ray  from  0  to  3  without  going  through  internal 
reflections  in  the  dielectric  shell.  Its  fle.ld  solution  is  given  in  (5.3) 
which,  of  course,  is  the  main  ccatrlbutlon.  The  second  type  contains 
rays  which  bounce  one  or  more  times  inside  the  shell  before  reaching  point  3. 
We  now  consider  Che  contribution  of  such  multiply  refracted  rays.  For  the 
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axial  Incldanca  caaa  (Flgura  23) ,  l«t  us  furthar  apc.clallca  taa  gaosatry 


T 


(5.20c) 


slab 


[a  -  ilz 
"  U  ♦  ij 


•xp(-j 2kQnb) 


Coaparing  (5. 20b)  vlt!'  (S.17),  wu  twtc  th«  efface,  of  the  dmltipl.y  refracted 
raya  is  accounted  for  by  replacing  T  by  When  the  condition  in 

(5. 20a)  la  not  aat,  ve  aust  avaluata  (5.19)  nuaarlcally.  Let  ua  define 
an  error  tera 


includlu'i  p  aultiply  refracted  raya 


»  lOOX.  (5,21) 


(1^0)1  including  all  aultiply  refracted  raya  j 
In  7igure  2?,  we  plot  2^  (including  no  aultiply  refracted  raya)  and 
(including  one  uuJtiply  lefracted  ray)  va.  u  for  b  ■  0,75  and 
(a  c)/b  «  A.  Several  obA'ervatlona  are  made,  (a)  For  coaaonly  uaed 
values  of  n  (between  1  and  3),  the  error  ia  13X  or  leaa  except  at 
resonancea.  A  r'^aonance  occurs  when  all  aultiply  refracted  raya  emerging 
from  the  alab  are  in  phaae  with  the  priaiary  ray,  Fcr  the  configu^raciou 
I?'  Figure  26.  the  resonance  condition  la 


(nk^b/iv)  -  (Inb/Xp)  ■  a  positive  integer  .  (5.22) 

For  an  obliquely  incident  ray  and/or  a  curved  slab,  Che  condit^cn  for 
resonance  is  rarely  satisfied.  Thus,  generally  speaking,  the  error  for 
neglecting  the  multiple  refraction  is  roughly  lOZ.  (b)  Errors  6^  aud  5^ 
have  about  the  same  order  of  aagnicude.  Thus,  the  iviclusion  of  the  first- 
order  aniltlply  refracted  ray  does  not  in  general  laprove  the  accuracy  of 
the  solution. 
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E.  Nu— rlcal  Reaulf 


(a)  Radowe.  Consldar  a  dlalaccrlc  apharlcal  radoaa  with  an  Innar 
radiua  20  a  unlfoni  thicknaaa  0^5  and  n  ■  /S,  aa  shown  in  Figura  28. 
Tha  y-polarisad  source  is  located  at  points  1,  3,  or  4,  and  ita  radiation 


field  in  the  B-plans  is  given  by 

"JV 

i^(r,e,4  ■  ir/2)  •  - — - —  [y  cos  (1.5  6)]  . 


(5.23) 


The  transaitted  field  is  calculated  froai  (5.3).  For  an  observation  point 
in  the  E-plane  and  in  the  far-field  sone  (x  ■  0  plana  and  c  in  Figure  21), 
we  may  rewrite  (5.3)  aa 

I^(r)  -  *-y-  [y  p^(e)l  ,  r^-  (5.24) 

where  r  Is  the  distance  from  point  0  to  point  3  (Figure  21).  Ue  plot  p^(e) 
as  a  function  of  d  in  Figure  28.  Generally  speaking,  the  radosie  iKKllfies 
the  radiation  field  gently,  as  expected. 

(b)  Lens  I  (Double  Concave).  Unlike  the  above  radome,  a  dielectric 
lens  nay  modify  the  incident  field  drastically.  Let  us  consider  Lens  I, 
drawn  approximately  to  scale  In  Figure  29.  The  source  is  2  X^  swey  from 
the  lens,  and  is  y-polarised.  Let  us  concentrate  on  the  field  in  the 
H-plane  (x-z  plane).  In  Figure  29,  %re  launch  4  rays  2*  apart.  The  outside 
ray  (at  9  ■■  6*)  suffers  total  reflection  at  the  second  face  of  the  lens, 
and  is  not  transnicted  into  the  free  space  region  (ve  Ignore  sniltiple 
refractions).  We  assume  the  Incident  field  from  the  source  is  confined 
to  a  cone  (a  beam) .  In  the  x-z  plane,  it  is  given  by 


if  e  <  9^ 

if  9  >  9^ 


(5.25) 
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0  ( DEG) 


.1  dielectric  lens:  Reo«etry 


In  th«  pr«s«nt  ca««  (Lent  I),  wt  choott  6^  >  5.5*.  After  trenntlttlon 
through  the  lent^  the  reyt  becoae  aore  divergent,  and  they  are  spread 
over  an  angular  region  of  about  6  <  53.5*.  which  it  contidcrably  wider 
than  the  incident  angular  spread  (6  <  5.5*).  At  an  obaarvatlon  point 
(r,0,4  ■  0),  we  axpreaa  the  tranaaitted  field  at 

I^(r.6.0)  -  (y  Q^r,e)l  .  (5.26) 

In  the  absence  of  the  lens,  ■  1  for  6  <  6^  and  ■  0  for  6  >  6^. 

t  3 

With  the  lens  present,  we  ploc  Q  (r,6)  as  a  function  of  6  for  r  *  2  x  10 

in  Figure  30.  Note  that  the  transaitted  field  is  auch  weaker  (12Z  or  less) 

than  the  incident  field  beer,  ise  of  the  wider  spread  of  the  transaitted  rays. 

For  the  present  case,  Q^(r,6)  is  only  very  weakly  dependent  on  r,  as 

long  as  r  >  100  X^.  Thus,  the  transaitted  field  in  (5.26)  in  the  far 

zone  is  approximately  a  spherical  wave  with  an  angulsr  pattevn  Q^.  Every 

transmitted  ray  has  two  foci.  Their  distances  behind  the  second  face  ^2 

of  the  lens  are  (>^21^  ^  ^*^22^  calculated  froa  (5.7). 

In  particular,  (<^21^  ^  normal  section  of  the  wavefront  in  the 

plane  of  incidence  (x-z  plane) ,  whereas  (•'22^”^ 

perpendicular  plane  (defined  by  the  y-axls  and  Che  ray  direction).  For 
the  third  ray  (incident  9  ■  4*)  in  Figure  29,  we  calculate  from  (5.7)  that 

(k^j^)”^  «  +0,21  Xq  ,  (k22>'^  •  +0.375  X^.  (5.27) 

Ue  mark  the  position  of  the  focus  corresponding  to  <21  by  a  cross  in 
Figure  29,  The  distance  between  A  and  the  cross  is  (<22)''*  trace 

of  the  two  sets  of  foci  is  shown  in  Figure  31.  They  are  curves  on  the 
two  caustic  surfaces  of  the  transmitted  wavefront  (intersection  of  caustic 
surfaces  and  Che  x-z  plane) . 
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Fi^iirt*  10.  far  ff<’M  n.itforn  riir(Hi|fli  lens  1  and  lens  II 


(c>  L«n»  II  (Doubl«  convx)  >  Th»*  gcoaatry  of  Lena  11  itnd  Its  ray 
plctura  ara  ahown  in  Flgura  32.  Tha  Incldant  fiald  ia  glvan  in  (S.25)  with 
6^  *  10*.  Tha  kfanacitted  fiald  in  rha  H-plana  ia  axpraaaad  in  (5.26),  whara 
Q^(t,0)  ia  again  vary  waakly  dapaudrnt  on  r  in  tha  far  tona  and  ia  plottad 
in  Figura  30.  Wa  nota  that  tha  incldant  10*-baaa  ia  nov  apraad  into  a 
67”*baaBi  after  tranaaiaaion  through  Lens  II.  Tha  paak  valua  of  tha 
tranaaittad  fiald  ia  about  24X  of  tha  incidant  field  at  tha  aaaa  far-field 
location.  Tha  cauatic  curvaa  ara  aho*m  in  Figura  31. 

(d)  Lana  III  ((kiovaxo-Concava) .  Tha  gaoaetry  of  Lena  III  ia  ahown 
in  Figura  33.  The  radii  <R^,R2)  of  the  lena  and  the  aource  distance 
satisfy  tha  lensMkar'a  aquation  in  (S.16),  so  that  the  tranaaittad  rays 
naar  the  axial  direction  ara  alaost  parallel  and  focua  at  a  point  at 
infinity  in  the  axial  direction.  The  incident  field  is  given  in  (5.25) 
with  0^  ■  12*.  Tha  H-plane  tranaaicced  field  expressed  in  (5.26), 

where  O^(r,0)  varies  drastically  froa  the  naar  field  aont  up  to  r  '  100 
as  seen  froa  Figure  34.  Beyond  r  '  100  X^,  the  beaa  becoaes  narrower,  and 
the  peak  becoaes  higher  with  the  increase  of  r.  It  is  wall-known  that 
the  exact  value  of  the  beam's  peak  (on  a  caustic  surface)  cannot  be 
predicted  by  the  present  geoaetrical  optics  theory.  It  can  be  calculated 
from,  for  e.g.,  the  Huygens-Green  formula  described  in  p.  107  of  [15]. 

The  caustic  curves  are  shown  in  Figure  35. 
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Figure  32.  Double-convex  spherical  dielectric  lens:  geo^try  and  ray  picture. 


e  (0E6J 

Figure  ?4.  H-plane  far  field  pattern  through  lena  III 


F.  Conclusions 


For  a  given  incident  field  in  (S.l),  a  geometrical  optics  solution  of 
the  transmitted  field  through  a  spherical  dielectric  shell  (Figure  22)  is 
given  in  (5.3)  through  (5.7).  The  present  solution  is  an  approximated  one. 
We  siunmarlze  its  validity  and  limitations  below 

(1)  Because  ray  techniques  are  being  employed  for  the  analysis 
presented  herein,  the  radii  of  curvature  should  be  large 
in  terms  of  the  incident  wavelength  in  order  for  the  results 
to  be  valid.  As  in  all  ray-optical  solutions,  our  solution 
fails  when  the  observation  point  is  near  the  caustic  surfaces 
of  the  transmitted  wavefront. 

(ii)  The  solution  given  in  (5.3)  and  (5.7)  is  valid  only  for 
spherical  shells,  and  when  four  points  (0,1, 2, 3,)  in 
Figure  22  are  coplanar.  These  two  limitations  can  be  relaxed. 

In  fact,  an  explicit  GO  solution,  (which  is  very  similar  to 
the  present  one,)  for  an  arbitrarily  curi/ed  shell,  and 
arbitrary  source  and  observation  points  is  given  in  [5]. 

(iii)  Cur  solution  includes  only  the  contribution  from  the  first- 
order  refracted  ray  (a  direct  ray  from  the  source  to  the 
obser'/atlon  point  withc'it  internal  refractions  inside  the 
dieiectri.c  shell).  It  has  been  estimated  that,  except  at 
"resonances"  which  rarely  occur  for  curved  shells,  the 
error  of  neglecting  higher-order  refracted  rays  is  roughly 
132  for  n  •  /c^  in  the  range  1  to  3.  It  should  be  pointed 
out  that  the  hlgher-ordei:  refracte-d  rays  can  be  calculated 
in  exactly  the  same  manner  as  the  first-order  refracted 
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ray  by  using  the  formulas  derived  in  this  work  and  in  [ 5 ] . 

It  is  only  a  matter  of  bookkeeping  and  computer  tino  in 
doing  the  numerical  calculations. 

(iv)  For  transmission  through  a  curved  dielectric  shell,  the 

major  field  contribution,  of  course,  comes  from  the  refracted 
rays  which  are  studied  In  this  paper.  However  other  diffraction 
processes  exist  wliich  may  not  be  conveniently  fitted  into  ray 
descriptions,  and  their  contribution  may  be  significant  under 
certain  conditions.  An  example  is  the  modal  fields  guided  by 
the  curved  dielectric  shell.  Efforts  to  understand  its  behavior 
in  canonical  problems  have  been  initlcted  [16]. 


VI.  ASSAY  COVESED  BY  SADQME 


In  chls  chapter,  ve  study  the  radiation  pattern  of  a  linear  array  of 
eleven  point  sources  enclosed  In  a  radone.  The  array  elenents  are  arranged 
symaetrlcally  along  the  x-axla  with  a  half-wavelength  spacing  as  shown  In 


Figure  36.  The  element  pattern  la  assumed  to  be  of  the  form 


- 


‘:5V 


[sin  ^  6  -(■  cos  ^  cos  6 


(6.1) 


To  calculate  the  array  pattern,  rays  are  traced  from  each  source  through 
the  radome  as  discussed  in  Chapter  I.  When  the  ray  tracing  Is  done,  there 
are  two  options  available  In  calculating  the  far  field:  (a)  we  can  calculate 
the  transmitted  field  over  a  planar  surface  just  outside  the  radome 
and  then  integrate  this  field  to  get  the  far-fleld  pattern;  (b)  we  can  push 
the  surface  to  Infinity  and  calculate  on  this  surface  directly,  thus 
avoiding  the  integration  step-this  Is  known  as  the  direct  ray  method. 

Though,  In  general,  both  methods  do  give  Identical  results,  Che  second 
method  becomes  less  accurate  If  the  far-fleld  point  happens  to  be  In  the 
vicinity  of  the  caustic  surfaces  of  the  transmitted  wavefront.  However, 
one  does  not  usually  confront  such  a  situation  In  practice.  The  details 
of  the  calculations  using  these  two  methods  may  be  found  in  [5]. 

Once  the  far-fleld  pattern  of  each  source  Is  obtained,  the  sum  pattern 
Is  obtained  by  adding  Che  patterns  of  the  individual  elements.  A  progressive 
phase  shift  Is  Introduced  among  Che  elements  to  enable  beam  scanning.  In 
Che  results  Co  follow,  we  present  Che  H-plane  (x-z  plane)  sum  patterns.  Thus, 
restricting  ourselves  to  this  plane,  the  Incident  field  In  (6.1)  reduces  to 
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ccj  6 


(6.2) 


Th«  total  transmitted  field,  auMsed  over  all  the  sources,  nay  then  be 
represented  by 


(6)  ,  e _ ^  Ie^  I  ^j(n-l)(ltdsine4a) 

sun  t  ^ ,  (6)  n 

n-4 


(6.3) 


where  |E  (0)|^  represents  the  magnitude  and  *^(6),  the  phase,  of  the  trans¬ 
mitted  field  due  to  the  source,  d  in  the  element  spacing,  and  a  Is 
the  Inter-element  phase  shift. 

Array  patterns  were  calculated  for  four  types  of  radomes  as  specified 
In  Table  I.  Two  of  the  radomes  were  paraboloids  and  two  were  tangent 
ogives . 

(^)  E*dome  A;  This  Is  a  paraboloidal  ladome  cf  relative  dielectric 
constant  2.S.  The  Inner  and  outer  surfaces  of  this  radome  were  generated 
by 

z/Xq  -  50  -  (x^  +  y^)/(8xj)  (6.4) 

and 


z/Xq  -  50.25 


(x^  +  y^)/(8.16xj)  , 


(6.5) 


respectively.  The  geometry  of  this  radome,  along  with  the  eleven  element 
array,  are  shown  In  Figure  36.  The  H-plsne  element  patterns  through  the 
radome  are  shown  In  Figure  37.  The  H-plane  swi  patterns  are  shown  In 
Figures  38-42  for  different  scan  angles  (the  scan  angles  Indicated  In  the 
figures  are.  the  values  for  the  beam  without  the  radome,  meaaured  from 
the  z-axls). 
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Figure  AO.  Same  as  Figure  38,  except  for  ^can  angle 


POLAR  ANGLE.  9  (DEGREES) 

Flgura  41.  S«at  Figure  38,  exc.«pc  for  acun  angle  •  -19.77*. 
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POLAR  ANGLE.  «  (DEGREES) 


Figure  42,  Same  ac  Figure  3S,  except  for  scan  angle 


SAHPE 

e 

r 

LENGTH 

(External) 

DIAMETER 

(External) 

Paraboloid 

2  5 

50.25  xj 

40.50  Xq 

A  •  j 

50.50  Xq 

40  Xq 

Tangent 

5.7 

14.6  Xq 

6.7  Xq 

Ogive 


WALL 

TRICmSS 


_ 

yz-o 

0.233  Xq 


Tlgur*  43  ahowi  tht  radoM  boras ight  arror  for  this  radoaa.  Tha  bora- 
sight  arror  is  tUa  difxaraoca  botwaan  tha  baas  saxisua  with  and  without 
tl.a  radosa  [  (baas  saxisua  with  radosa)  -  (baas  saxisus  without  radosa)  ] . 

It  say  alao  ba  obsarvad  fros  Figuraa  38-42  that,  as  tha  scan  angle 
incraaaas,  tha  gain  dacraaaas  bacauaa  of  tha  dacraaaa  in  tha  affactiva 
apertura  with  incraaaing  scan  angla.  Figura  44  shows  this  gain  variation. 
Obsarva  that  in  tha  prasanca  of  tha  radosa  tha  saxis'js  fiald  strangth  drops 
such  faster.  Also,  around  the  tip  region  where  tha  curvature  is  tha  highest, 
the  saxisia  field  strength  decreases  faster,  and  than  lavals  off  at  higher 
sran  angles. 

(b)  Redone  B;  This  is  also  a  paraboloidal  radons  of  ■  2.S,  but 

the  outer  surface  in  this  radone  was  configured  so  as  to  nnke  tha  curvatvra 
in  the  tip  region  very  close  to  the  critical  curvature  (sea  Figure  30,  [5]), 
thus  resulting  in  a  greater  focussing  action  in  tha  s-diraction.  Tha  inner 
surface  was  obtained  from  (6.4)  and  the  outer  surface  was  obtained  fron 

s/Xq  «  50. 5C  -  (x^  +  y^)/(7,925XQ)  .  (6.6) 

The  geonetry  of  the  radoce  is  shown  in  Figure  45.  Obsarva  that  tha  thickness 
is  naxlnun  along  the  t-fiXis  and  reduces  gradually  to  saro  thickness  at  tha 

base. 

The  H-plane  sun  patterns  are  shown  in  Figures  46-49,  and  the  boresight 
error  curve  Is  presented  In  Figure  50.  The  focussing  action  of  this 
radone  may  be  clearly  observed  in  Figure  46  -  the  axial  field  strength  is 
1.4  tines  as  large  as  that  without  tha  radone. 


SCAN  AN<iLE  ( DEGREES ) 


Figur*  43.  Bor«3lght  arror  for  radooo  A, 


•i 


POLAR  ANGLE,  9  (DEGREES) 


Figure  46.  sua  pactern  through  radoM  B.  Boaa  at  0*. 


80RES1GHT  ERROR  (DEGREES) 


(c)  C  ^od  Dt  Th«t«  ar*  t«nt«nt  o|lv«  radoMi  in  ■Inaila*. 

Th«  ••OMtrlcal  conf Iguration  la  givan  la  flgura  SI.  fUdoM  C  haa  •  S.7 
atul  radoM  ^  ■  9.3.  It  may  ba  obaarvad  fro*  Flgura  SI  that  tWa  radcM 

haa  a  dlacontlnulty  In  tha  a-dlractlon.  In  pracclca,  thla  dlacontiouona  tip 
raglon  la  raplacad  by  a  a*ooth  cunrad  aurfaca.  Hotfavar,  In  our  eaicuiatlona, 
chit  nodlficatlon  uaa  not  Introducad.  lha  tlald  along  tha  a-axla  vaa 
calculatad  by  Intarpolatlon. 

Tho  K^plana  SM  pattarna  for  radoma  C  ara  praaantad  In  Plguraa  S2-S7, 
and  tha  corvaapondlng  boraalght  atror  curv-a  la  ahoam  In  Flgura  i»3.  Ilta 
pattama  for  radooa  D  ara  in  Flgurca  S9-64.  The  Soraalght  error  curve  la 
ahown  In  Flgura  6S.  It  nay  i»a  obaarvad  fron  thaaa  flguraa  that  tha  pattama 
for  anall  acan  anglaa  ara  not  i4:anarally  good.  A.lac,  tha  boraalght  atror 
at  anall  acan  anglaa  xa  quite  high.  Thla  behavior  la  not  uaually  obaarvad 
In  practical  radooMa  for  tvo  raaaona: 

(li  In  practice,  the  ancannx  «  uaually  r  reflactor  ontanna  or  a 
planar  array  of  coniidarably  large  nunbar  of  elaMnta»  uharaaa 
In  our  nodal,  va  have  only  alavtu  point  auurcaa  Xq/2  apart. 

If  one  Includaa  aufflclant  nuaboi  of  point  aourcaa  ao  aa  to 
faithfully  roproaant  a  given  practical  antenna,  there  vlll  be 
conalclarabie  rhengax  in  c^a  aua  pactema  and  Che  boreelght  error 
curvea. 

(11)  Tha  eacond  raaaon  haa  aooathlug  to  do  vlth  tha  tip  raglon. 

Figure  64  ■!iOwo  a  ruigant  ogive  end  »  paraboloidal  radone  of 
tha  aana  alia.  It  nay  be  obaerved  that  the  tip  region  In  tho 
rangact  ogive  la  vary  ataap  wharaaa  In  tha  paraboloid  tha 
aurfaca  changes  anra  anoothly. 

In  Qaniral,  tha  pat  cams  Inprova  at  largar  acan  anglaa,  aa  doea  the 


boraalght  error. 
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SUM 


POLAR  ANGLE.  9  (DEGREES) 

Flgura  S7.  Saa*  Figure  52,  except  for  scan  angle  ■  -27.64*. 
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SORESIOHT  ERROR  (OEOREES) 


-30  -20  -10 


0 


10  20  30 


POLAR  ANGLE  .  B  (DEGREES) 

Flfurt  59.  H-pl«o«  pattern  through  radoM  D.  Baaa  at  0*. 
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SCAN  ANGLE  (DEGREES) 


Flgurv  65.  Bor««lghc  error  for  raUoM  D. 
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APPENDIX  A 


DERIVATION  OF  EQUATION  (5.6) 


Referring  to  Figure  21,  the  curvature  natrix  fA],  [10] of  the  incident 
spherical  wave  may  be  expressed  as 

0  ■ 

,  .  (A.l) 

)  a"\ 

The  curvature  matrix  of  the  Inner  surface  of  the  radome  (I^)  Is 


,-l 


(A.  2) 


From  the  knowledge  of  these  two  curvature  matrices,  the  curvature  matrix  of 
the  refracted  ray  1-2  nay  be  expressed  as  [A] 


qJ  -  (0j)"^[6jQjeJ/n  +  (cos  oj  -  cos  aJ/n)Q 
where  6^  and  6^  are  coordinate  transformation  matrices  given  by 


(A.  3) 


cos  a. 


(A.  A) 


cos  u. 


(A.  5) 


Simplifying  (A. 3),  we  obtain 


(A.  6) 

wh«r*  and  k^2  defined  in  (5.7). 

The  curvature  matrix  of  the  wavefront  Incident  at  point  2  it  given  by 

Qj  -  +  b  !l‘^  (A. 7) 

where  T  is  the  identity  matrix.  The  curvature  matrix  of  the  transmitted 
wavefront  at  2  can  be  calculated  similarly  to  that  at  1.  The  final  result 
is 


where  *^22  defined  in  (5.7). 
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Coluidiua,  OH  43212 

Univaraity  of  Colorado 
Dapartaent  of  Elactrical  Englnaaring 
ATW:  Prof.  L.  Lawin 
Bouldar,  CO  80302 

Ganaral  Dynaaica  Corporation 
Elactronica  Dlviaion 
ATTM:  Dr.  G.  Trlcolaa 
P.  0.  Box  81127 
San  Diago,  CA  92138 

Coaaandar 

Air  Forca  Avionica  Laboratory/IWf 
ATTN:  Mr.  Harold  Vabar 
Mr.  Alan  Bluaa 

Vright-Pattaraon  Air  Forca  Baaa,  OH  45433 

a 

Offica  of  Naval  Haaaarch 
ATTN:  Dr.  H.  Nullanay/Coda  427 
800  North  Quincy  Straat 
Arlington,  VA  22217 

R.  C.  Hanaan,  Inc. 

Box  215 

Taraana,  CA  91356 
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Harry  Diaaond  Laboratorlaa 

ATTN:  Dr.  Hovard  S.  Jonaa.  Jr./AHXDO-RAS 

2800  Powdar  Mill  Itoad 

Adalphl,  MD  20783 

Mlaalla  Intalliganca  Agaecy 

US  Army  Mlaalla  CoMand 

ATTN:  Mr.  R.  Thoapaon  (MSMI/TPE) 

Radatona  Araanal,  AL  35809 

Rughaa  Aircraft  Coapany 

Radar  Syataaa  Group 

ATTM:  Dr.  V.  H.  Kuaaar/MSR2-Al02 

P.  0.  Box  92426 

Loa  Angalaa,  CA  90009 

Mr.  Gordon  Lovlcc 

Offlca  of  Naval  Raaaarch  Branch  Office 
536  S.  Clark 
Chicago.  XL  60605 

Prof.  S.  W.  Laa 
Electroaagnaclca  Laboratory 
Dapartaant  of  Elactrlcal  Engineering 
Unlvaralty  of  Illlnola 
1406  W.  Graan  St. 

Urbana,  Illlnola  61801 

Prof.  R.  Mlttra 
Elactroaagnatlca  Laboratory 
Dapartaant  of  Electrical  Engineering 
Unlveralty  of  Illlnola 
1406  V.  Green  St. 

Urbana.  Illlnola  61801 

Mr.  H.  Lawler 

Dapartaant  of  Elactrlcal  Englnaarlng 
Unlvaralty  of  Illlnola 
1406  V.  Graan  St. 

Urbana.  Illlnola  61801 
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